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Abstract 

We present a study of r~ — y -k~ K®K^(-k )v t and r~ — > K~ K%K < i{jv°)v r decays using a dataset of 
430 million r lepton pairs, corresponding to an integrated luminosity of 468 fb _1 , collected with 
the BABAR detector at the PEP-II asymmetric energy e + e~ storage rings. We measure branching 
fractions of (2.31 ± 0.04 ± 0.08) x 10~ 4 and (1.60 ± 0.20 ± 0.22) x 10~ 5 for the r~ -»■ ty~ K° s K° s v T 
and t~ — ¥ tt~~ K$ K®n° v T decays, respectively. We find no evidence for r~ — \ K~ K%K%v T and 
r~ — > K~ K%K®n°v T decays and place upper limits on the branching fractions of 6.3 x 10 -7 and 
4.0 x 10~ 7 at the 90% confidence level. 

PACS numbers: 13.35.Dx, 14.60.Fg 



The t lepton can be used as a high-precision probe of 
the Standard Model (SM) and models of new physics. A 
recent BABAR paper, for example, presented a search for 
CP violation by measuring the decay-rate asymmetry of 
t~ — s- ir~K®v T decays [ij. One of the backgrounds in 
that analysis is t~ — > 7r~ K^K^i> t , which has a large un- 
certainty in the branching fraction Q. The uncertainty 
in the background from r~ — > n~ K®K®i/ T decays was not 
a limitation of the decay-rate asymmetry measurement, 
but an improved measurement of the branching fraction 
and an understanding of the decay dynamics will be re- 
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quired for a future measurement at a high-luminosity B- 
factory. 

This paper presents measurements of the branching 
fractions of t~ — > 7r~ K^K^tt )^ decays and the first 
search for t~ — > K~ K^K^tt )^ decays. In this work we 
use the K®— > 7r + 7r~ decay mode. Here and throughout 
the paper, charge conjugation is implied. 

Previously, ALEPH and CLEO measured the r~ — > 
tt~ K^K^ist branching fraction to be (2.6 ± 1.0 ± 0.5) x 
10~ 4 3] and (2.3 ± 0.5 ± 0.3) x 10" 4 0], respectively. 
ALEPH set an upper limit on the — > ir~ K®K®ir Q v T 
branching fraction of 2 x 10~ 4 at the 95% confidence level 

1- 

The present analysis uses data recorded by the BABAR 
detector at the PEP-II asymmetric-energy e + e~ collider, 
operated at center-of-mass (CM) energies of 10.58 GeV 
and 10.54 GeV at the SLAC National Accelerator Labora- 
tory. The BABAR detector is described in detail in Ref . Q ■ 
In particular, charged particle momenta are measured 
with a five-layer double-sided silicon vertex tracker and 
a 40-layer drift chamber, both within a 1.5 T supercon- 
ducting solenoidal magnet. Charged kaons and pions are 
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separated by ionization (dE/dx) measurements in the sil- 
icon vertex detector and the drift chamber in combina- 
tion with an internally reflecting Cherenkov detector. An 
electromagnetic calorimeter made of thallium-doped ce- 
sium iodide crystals provides energy measurements for 
electrons and photons, and an instrumented flux return 
detector identifies muons. Based on an integrated lu- 
minosity of 468 fb -1 , the data sample contains approxi- 
mately 430 million r-pair events. 

Simulated event samples are used to estimate the se- 
lection efficiency and purity of the data sample. The 
production of r pairs is simulated with the KK2F Monte 
Carlo (MC) event generator @. Subsequent decays of 
the r lepton, continuum qq events (where q = u,d,s,c), 
and final-state radiative effects are modeled with Tauola 
and EvtGen @, JETSET @, and PHOTOS Q3, re- 
spectively. Passage of the particles through the detector 
is simulated by Geant4 [111 ]. 

The t~ — >• 7r~ K®K2,v T decay is simulated with Tauola 
using t~ -> K*~K v T . The t~ ->■ 7r~ K°K°n°v T de- 
cay is simulated with EvtGen using t~ — > K*~ K tt°v t 
and r~ — > K*°K a Tr~v r . As we later show, the t~ — » 
K*~K°u T and t~ -> K*-K°tt°i> t have a if* (892) me- 
son that is observed in the ir~ K® channel, and the 
t~ — s- K*°K°ir~v T has a if* (892) meson that is observed 
in Tr°K® channel. 

The r pair is produced back-to-back in the e + e~ CM 
frame. As a result, the decay products of the two r lep- 
tons can be separated from each other by dividing the 
event into two hemispheres - referred later as the "sig- 
nal" hemisphere and the "tag" hemisphere - using the 
plane perpendicular to the event thrust axis jl~2j]. The 
event thrust axis is calculated using all charged particles 
and all photon candidates in the entire event. 

We select events with one prompt track and two K® 
— > tt + 7r~ candidates reconstructed in the signal hemi- 
sphere, and exactly one oppositely charged prompt track 
in the tag hemisphere. A prompt track is defined to be 
a track with its point of closest approach to the beam 
spot being less than 1.5 cm in the plane transverse to the 
e~ beam axis and less than 2.5 cm in the direction of the 
e~ beam axis. Tracks consistent with coming from a K® 
or A decay, or from a 7 conversion are not considered to 
be prompt tracks. The components of momentum trans- 
verse to the e~ beam axis for each of these two prompt 
tracks must be greater than 0.1 GeV/c in the laboratory 
frame. A K® candidate is defined as a pair of oppositely 
charged pion candidates with an invariant mass between 
0.475 and 0.525 GeV/c 2 (see Fig. [J); furthermore, the dis- 
tance between the beam spot and the n + n~ vertex must 
be at least three times its uncertainty (the di-pion pair 
will be referred to as the "K® candidate daughters"). 

The charged hadron must be identified as a charged 
pion or a charged kaon. The efficiency for selecting 
charged pions and kaons is approximately 95% and 90%, 
respectively. The probability of mis-identifying a charged 
pion (kaon) as a charged kaon (pion) is estimated to be 
1% (5%). 



The charged pion and kaon samples are divided into 
samples with zero and one 7r° mesons. Events with two or 
more ir° mesons are rejected. The 7r° candidate is recon- 
structed from two clusters of energy deposits in the elec- 
tromagnetic calorimeter that have no associated tracks. 
The energy of each cluster is required to be greater 
than 30MeV in the laboratory frame, and the invariant 
mass of the two clusters must be between 0.115 GeV/c 2 
and 0.150 GeV/c 2 . The clusters in the electromagnetic 
calorimeter that are not associated with a ir° candidate 
are ignored in the analysis. 

To reduce backgrounds from non- r-pair events, we re- 
quire that the momentum of the charged particle in the 
tag hemisphere is less than 4 GeV/c in the CM frame 
and be identified as either an electron or a muon. For 
momenta above lGeV/c in the laboratory frame, elec- 
trons and muons are identified with efficiencies of ap- 
proximately 92% and 70%, respectively We also 
require the magnitude of the event thrust to be between 
0.90 and 0.995. 

The invariant mass of the charged hadron and the two 
K® mesons is required to be less than 1.8 GeV/c 2 . For 
t~ — s- 7r~ K^KgTT ^ decays, we do not include the n° 
in the mass calculation. The ii~K®K2. invariant mass 
is shown in Figs. [5] and [3] The tt~ K®K®ir° invariant 
mass is also shown in Fig. [3] We also require that 
the pseudomass to be less than 1.9 and 2.1 GeV/c 2 for 
the t~ — > ir~ K°KgU T and r~ -» n~ K®K®-k°v t sam- 
ples, respectively (the 7r° meson is included in the pseu- 
domass calculation). The pseudomass is defined to be 
A/pscudo = yjMl + 2(v/i - E h ){E h - P h ) where E h and 
Ph are the energy and magnitude of the momentum of 
the hadronic final state in the laboratory frame [l4| . 

The invariant mass distribution predicted by the MC 
for the hadronic final state particles and for their combi- 
nations do not perfectly describe the data. In particular, 
the peak of the (ir~~ K®K®) invariant mass distribution in 
the MC is found to peak approximately 5% lower than 
the peak observed in the data. To improve the modeling 
of the data we have weighted the t~ — > tt~ K®K®v r in 
Tauola using the Dalitz plot distribution for the tt~K® 
invariant mass (shown for the data sample in Fig. [2]). The 
weighting function is from a two-dimensional (9 x 9) ma- 
trix using M 2 (ir~ K®) with both ir~K® combinations (the 
matrix is constructed to be symmetric). The weighted 
events are used in all the mass plots and we observe an 
improvement in the modeling of the data. 

The branching fractions of the two charged pion modes 
are determined simultaneously to take into account the 
cross feed of each decay mode into the other sample. The 
branching fraction is 

Bj =E l ^\N^-N^ d )/(2N TT ) 

where j represents the t~ — > tt~ K®K®iy T and t~ — > 
tt~ K^K^tt ^ decay modes; i represents the (n~KgKg) 
or (tt-X^tt ) reconstruction modes; N? ata and A^ kgd 
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are the number of data and background events in the i-th 
data sample; e _1 is the inverse of the selection efficiency 
matrix (ey is the probability to select an event of type j 
with the selection criteria i); and N TT is the number of 
r-pair candidates. 

The columns in Table U give the number of data and 
background events for each reconstruction mode. Ta- 
ble Q] also gives the selection efficiency matrix, where the 
horizontal row gives the efficiency for selecting the true 
decay for each reconstructed mode. For example, the ef- 
ficiency for selecting a true r~ — > 7T~KgK®v T decay is 
(4.93 ± 0.02)% and (0.21 ± 0.01)% with the (tt- K°K° S ) 
and (tt~ K^K^tt ) selection criteria, respectively. 

We measure the t~ — > n~K®K®v T and r~ — > 
7T~ K®K®ir v T branching fractions to be 



B(j- - 



-> n-K°K°u T ) = (2.31 ± 0.04 ± 0.08) x 10~ 4 
n~K°K°ir v T ) = (1.60 ± 0.20 ± 0.22) x 10~ 5 
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FIG. 1: The invariant mass of the two — > it + tt~ candi- 
dates in the r~ — > -k~ K®K®v T (top) and r~ —¥ n~ K®K®-k°v t 
(bottom) samples after all selection criteria have been ap- 
plied. The points are data and the histograms are the 
prediction of the Monte Carlo simulation. For both plots, 
the white histogram represents r~ — >• K*~ K°v T decays, the 
blue and beige histogram shows the r~ — > K*~ K and 
T - _^ K*°K°ti-v t (r~ ->• ■K-K%K 3 n°v T ) decays, respec- 
tively. The red histogram is the qq background. 



where the first error is statistical and the second is sys- 
tematic. The statistical correlation parameter for the two 
measurements is found to be —0.21. The results have 
been corrected for the K® — > ir + ir~ branching fraction 

0. 

The systematic uncertainties (see Table A} are divided 
into the selection efficiency, background, and common 
systematic components. The uncertainties on the ele- 
ments of the efficiency matrix only include the errors spe- 
cific to that decay and selection criteria. Uncertainties 
that are common to all matrix elements are included in 
the common systematic errors. 

The efficiency for selecting t~ — > n~ K a s K®v T events is 
found to be (4.93±0.02)% and (0.21±0.01)% for the sam- 
ples with zero and one ir° candidate, respectively. The 
uncertainty on the first efficiency is from the MC statis- 
tical error. The uncertainty on the second efficiency also 
includes the MC statistical error and an error that takes 
into account the uncertainty for finding a fake ir° meson 
in t~ —> ir~ K®K®v T decays. The uncertainty for find- 
ing a fake ir° is estimated to be 6% and is determined 
by comparing the number of t~ — > tt~ K®K®v T decays 
that have two neutral clusters in the data and MC sam- 
ples where the invariant mass of the two neutral clusters 
must not be near the 7r° mass. 

The efficiency for selecting r~ — > tt~ K^K^ir ^ events 
is found to be (3.04 ± 0.10)% and (2.65 ± 0.09)% for the 
samples with zero and one ir° candidate, respectively. 
The uncertainties include the MC statistical error and 
an uncertainty for the 7r° identification. The uncertainty 
for identifying a 7r° meson is estimated to be 3% based 
on studies involving a variety of data and MC control 
samples. We observe that the efficiency for selecting 
t~ — > 7r~ K a s K®ir a v T decays with and without a ir° is 
approximately equal, and hence we assign a 3% uncer- 
tainty on the efficiency for selecting t~ — > ir~ K2,K2.tt v T 
decays without reconstructing the 7r° meson. 

The background in the charged pion modes is predicted 
by the MC simulation to be entirely from e + e~ — > qq 
events. The background in the charged kaon modes is 
cross-feed from the charged pion modes where a charged 
pion is mis-identified as the charged kaon. The back- 
ground in the charged pion sample is confirmed with 
data and MC simulation control samples. The control 
samples are created using the nominal selection criteria 
except that the invariant mass and pseudomass require- 
ments are reversed to eliminate the r-pair events and 
enhance qq events. The ratio of selected events in the 
data to MC control samples is found to be consistent 
with unity within 15% for both r~ — > ir~ K° s K° s v T and 
t~ — » 7r~ KsK®Tr°v T samples. The 15% value is added 
to the MC statistical uncertainty of the number of back- 
ground events. 

A number of systematic uncertainties are common to 
both the r" -> ■k~K°K°u t and r" -> n~K°K°w°v T 
branching fractions measurements. They can be catego- 
rized into two components: tracking and particle identi- 
fication reconstruction uncertainties, and topological se- 
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FIG. 2: The Dalitz plot of the (tt~ K°) system, and the (tt K°K^), (it K°) and (ifsiQ) invariant mass distributions for 
events that pass the r~ — > 7T~K° K%v T selection criteria. There are two entries per event in the Dalitz plot. The points are 
data and the histograms are the prediction of the Monte Carlo simulation. The signal decays are represented by the white 
histogram (r~ -> K*~ K°u T ). The beige histogram shows the r~ -> K*~ K°ir°i' r and r~ -> K*°K°tv-v t (t~ -» tt~ K s K° s tt°is t ) 
decays, The red histogram is the qq background. The mass plots use r~ — > n~ K®K®v T events that have been weighted based 
on the Dalitz plot distributions in the top left plot. 



lection uncertainties. 

The tracking and particle identification reconstruction 
uncertainties include the uncertainty on the track recon- 
struction efficiency (0.5%). They also include the un- 
certainties on the efficiencies of the particle identifica- 
tion algorithms: lepton identification (combined electron 
and muon) (1.6%), charged pion particle identification 
(0.5%), and K° s identification (1.8% for two K°). The 
particle identification algorithms used in this work are 
based on standard BABAR routines and the uncertain- 
ties are determined using control data and MC samples 
5, 15]. The uncertainty on the efficiency for selecting 
7T° mesons is included in the elements of the selection 
efficiency matrix. 

The topological selection uncertainties include a 2% 
uncertainty associated with the topological selection cri- 
teria that impose requirements that the prompt tracks be 
associated the primary vertex. Also included is the un- 
certainty in the product of the luminosity multiplied by 
the e + e~ — s- t + t~ cross section (1%). The weighting of 
the invariant mass distribution for the t~ — s- it~ K®K®v t 



MC decays does not change the number of events and, 
hence, does not alter the measured branching fractions. 

In Fig.dwe plot the (^K° S K°), (n"K°), and (K°K°) 
invariant mass distributions. The contribution of the 
if* (892) resonance (K* — > tt~ K°) is observed in the 
(n~Kg) invariant mass plot and the Dalitz plot in Fig. [2] 

The r~ — » 7r~ K a s K a s v T branching fraction is in good 
agreement with the previous measurements of (2.6±1.0± 
0.5) x 10~ 4 [| and (2.3 ± 0.5 ± 0.3) x 10~ 4 4j. The the- 
oretical prediction for the r~ — > tt~ K2,K®v t branching 
fraction is 4.8 x 10 -4 [lq ]. Decays involving a pion and 
two kaon mesons can have contributions from both ax- 
ial and vector currents at the same time, and the vec- 
tor contribution for r~ — s- tt - K®K®v t is estimated to be 
1.4 x 10~ 4 . 

Assuming isospin symmetry (l7j and using other mea- 
surements, we can estimate the r~ — > n~KgK^u r 
branching fraction. The t" — > n~ K° K°v T and t~ — > 
tt~ K + K~ v T branching fractions are equal if isospin is 
an exact symmetry (the t~ — > iv~KgKgV T and t~ — > 
tk~K\K\v t branching fractions are also equal). Hence 
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TABLE I: Results for the charged pion decays. The background events are primarily qq events. 



Decay mode 


T — > TT K S K S U T 


T — >• 7T K S K S 1\ ■ V T 


Branching fraction 


(2.3f ± 0.04 ± 0.08) x 10~ 4 


(1 .60 ± 0.20 ± 0.22) x 10~ 5 


Events 






Data 


4985 


409 


Estimated background 


98 ± 17 


35 ± 7 


Selection efficiency 








(4.93 ± 0.02)% 


(0.2f ± 0.01)% 


— — T^O T^O 

t — y 7r KsK s icv T 


(3.04 ± 0.f0)% 


(2.65 ± 0.09)% 


Fractional systematic errors 






Selection efficiency 


0.008 


0.12 


Background 


0.004 


0.04 


Common systematics 


0.034 


0.03 


Total 


0.035 


0.13 



B(t~ ir-K^Klvr) = B{t- ir-K+K-v T ) - 
2B(t~ — > -k~ K®K®v T ) and we obtain 

B(r~ -> n-K° g K a L v T ) = (9.8 ± 0.5) x 1CT 4 

where B{r~ -> -k-K+K~v t ) = (14.4 ± 0.4) x lO" 4 ]! 
based on measurements from BABAR [18| and Belle (l9j . 
The prediction is in good agreement with the branching 
fraction measured by ALEPH of (10.1 ± 2.3 ± 1.3) x 10~ 4 
i 

The MC simulation uses t~ -> K*~ K°tt®v t and t~ 
K*°K°tt~ ~v T decays. A model based on a phase space dis- 
tribution of the final state particles does not describe the 
(ir~ K®) invariant mass distribution. The relative contri- 
bution ofr~^ K*~K°tt°v t to r"-> K*°K°ir~v T decays 
is determined to be (0.17 ± 0.03) by simultaneously fit- 
ting the (ir~K®) and (tt K®) invariant mass distributions 
(see Fig. [3J. The predicted Monte Carlo distributions 
are fit to the data spectra after the subtraction of the 
t~ — > 7r~ K®K®is T and background events. The normal- 
izations of the two modes are varied with the constraint 
that the values be positive numbers. If we do not in- 
clude the r~ — > K*®K®tt~v t decay, then we observe a 
disagreement between the data and MC samples in the 
lower-mass and higher-mass regions of the M{it~K2.) and 
M(ir°Ks) distributions, respectively, in Fig. [3j 

The same criteria are used to select t~ — > K~ K a s K a s v T 
and r~ — > K~ K®K2.n°v T decays except that the charged 
track is required to be a kaon. The numbers of events 
are given in Table [TT] and found to be consistent with 
the estimated background prediction. The background is 
almost entirely due to cross feed of t~ — > tt~ K®K®is T and 
t~ — > 7r~ K®K®ir°i/ T decays and very little background 
from qq events. The branching fractions are determined 
for each channel independently and used to place upper 



limits on the branching fractions of 

B(t~ -> K-K Q s K° s v T ) < 6.3 x 10~ 7 

B(t~ -)• K~ K s K° s -k°v t ) < 4.0 x 10~ 7 

at the 90% confidence level. 

The t~ K-K°K°v T and r~ K~ K + K~ v T 
branching fractions are also predicted to be equal assum- 
ing isospin symmetry. The t~ — > K~ K + K~ v T branch- 
ing fraction is (2.1 ± 0.8) x 10 -5 based on measure- 
ments from BABAR [HI and Belle 0. BaBAR finds that 
a t~ —> K~(j)v T contribution can account for all of the 
t~ — s- K~K + K~v T decays. This suggests that the t~ — > 
K-K°K°v T and consequently, the t~ ^ K~K°K °ir°u T 
branching fractions should be small in the limit of isospin 
symmetry. 

In summary, we have measured the branching frac- 
tions of the r~ -» n~K®K®i/ r and t~ -> ir~ K®K°ir°v T 
decays to be (2.31 ± 0.04 ± 0.08) x 10~ 4 and (1.60 ± 
0.20 ± 0.22) x 10~ 5 , respectively. The r" -> n~K°K°u T 
and t~ — > ir~ K®Kgir°i/ T decays can be modeled using 
r~ -> K*-R Q v T , and r~ K*-K°it°v T and r~ 
jf*0jf0 7r - I y T decays, respectively. The r~ -> ^"ifgi^g^r 
branching fraction is a significant improvement on the 
previous measurements and the t~ — > 7r~-R"gi ; d?7r i/ T 
branching fraction is the first measurement. In addition, 
we place the first upper limits on the branching fractions 
of6.3xl0~ 7 and4.0xl0" 7 on the t~ -t K~K ° s K° s u T and 
t~ — » K~ K a s K° s ir°v T decay modes at the 90% confidence 
level. 
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TABLE II: Results for the charged kaon decays. The background events are primarily qq events. 



Decay mode 


r"-> K -R% K%v T 




Branching fraction 


(1.9 ±3.0 ±0.3) x 10" 7 


(1.5 ± 1.8 ±0.1) x 10" 7 


Limit (90% C.L.) 


6.3 x 10" 7 


4.0 x 10" 7 


Events 






Data 


23 


1 


Estimated background 


20.0 ±0.5 


0.15 ±0.02 


Selection efficiency 


(3.85 ± 0.04)% 


(1.37 ±0.03)% 
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